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Abstract. Spectrum of the neutrinos produced in the massive photon and longitudal plasmon decay process 
has been computed with four levels of approximation for the dispersion relations. Some analytical formulae 
in limiting cases are derived. Interesting conclusions related to previous calculations of the energy loss in 
stars are presented. High energy tail of the neutrino spectrum is shown to be proportional to exp(-E/kT), 
where E is the neutrino energy and kT is the temperature of the plasma. 
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Thermal neutrino loses from plasma are very important 
for stellar astrophysics [TJd] ■ Plasmon decay is one of the 
three main reactions. Extensive calculations for these pro- 
cesses were done by group of Itoh [3|l4 .5.6.7,8.9, lOtfTT] . 
Other influential article include P^ITWilfTMTIllTflfTSlfTgi 
[20] . Meanwhile, our abilities to detect neutrinos has grown 
by many orders of magnitude, beginning with 1.4 tonne 
experiment of Reines&Cowan up to the biggest exist- 
ing now 50 kt Super-Kamiokande detector [22]. Recently, 
"GADZOOKS!" upgrade to Super-Kamiokande proposed 
by Beacom&Vagins [23] attract attention of both experi- 
mental and theoretical physicists. At last one new source 
of the astrophysical antineutrinos is guaranteed with this 
upgrade, namely Diffuse Supernova Neutrino Background 
^ ' [2%1[24]. Pre-supernova stars will be available to observa- 
' tions out to ^2 kiloparsecs [5S]. This technique is the 
only extensible to megaton scale [25]. Memphys, Hyper- 
Kamiokande and UNO (Mt-scale water Cherenkov detec- 
tors cf. e.g. [55]) proposals now seriously consider to add 
GdCl3 to the one of the tanks with typical three-tank 
design [27]. Recently, the discussion on the geoneutrino 
detection [28j . increased attention to the deep underwa- 
ter neutrino observatories [29] with target mass 5-10 Mt 
[25] and even bigger [30]. It seems that (anti)neutrino as- 
tronomy is on our doorstep, but numerous astrophysical 
sources of the v's still are not analyzed from the detection 
point of view. 

Detection of the solar 31,32,33,34] and supernova neu- 
trinos [35,36,37,38] was accompanied and followed with 
extensive set of detailed calculations (see e.g. [31^l4T)ir4"Tl 
42,43,44 and references therein as a representatives of 
this broad subject) of the neutrino spectrum. On the con- 
trary, very little is known about spectral neutrino emission 
from other astrophysical objects. Usually, some analytical 
representation of the spectrum is used, based on earlier 
experience and numerical simulations, cf. e.g. [45] . Whilc 
this approach is justified for supernovae, where neutrinos 



are trapped, other astrophysical objects are transparent 
to neutrinos, and spectrum can be computed with an ar- 
bitrary precision. Our goal is to compute neutrino spec- 
tra as exact as possible and fill this gap. Plasmaneutrino 
process dominates dense, degenerate objects like red giant 
cores 46 , cooling white dwarfs 47J including la supernova 
progenitors before so-called ,,smoldering" phase [15] ■ It is 
also important secondary cooling process in e.g. neutron 
star crusts [49] and massive stars [50] . Unfortunately, ther- 
mal neutrino loses usually are calculated using mcthods 
completely erasing almost any information related to the 
neutrino energy £ v and directionality as well. This infor- 
mation is not required to compute total energy Q radiated 
as neutrinos per unit volume and time. From experimen- 
tal point of view, however, it is extremely important if 
given amount of energy is radiated as e.g. numerous keV 
neutrinos or one 10 MeV neutrino. In the first case we 
are unable to detect (using available techniques) any tran- 
sient neutrino source regardless of the total luminosity and 
proximity of the object. In the second case we can detect 
astrophysical neutrino sources if they are strong and not 
too far away using advanced detector which is big enough. 

Few of the research articles in this area attempt to es- 
timate average neutrino energy [16,17,52,53 computing 
additionally reaction rate R. Strangely, they presented fig- 
ures and formulae for Q/R instead of ^Q/R. This gives 
false picture of real situation, as former expression gives 
{£ v + £p). Obviously, we detect neutrinos not v-v pairs. 
hQ/R do not give average neutrino energy, as in general 
neutrino and antineutrino spectra are different. As we will 
see only for longitudal plasmon decay neutrinos energies 
of neutrinos and antineutrinos are equal. However, dif- 
ference in all situations where thermal neutrino loses are 
important is numerically small and formula: 



(1) 



is still a "working" estimate. 
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Mean neutrino energy is useful in the purpose of quali- 
tative discussion of the detection prospects/methods. Quan- 
titative discussion require knowledge of spectrum shape 
(differential emissivity dR/d£ v ). High energy tail is par- 
ticularly important from an experimental detection point 
of view. Detection of the lowest energy neutrinos is ex- 
tremely challenging due to numerous background signal 
noise sources e.g. 14 C decay for £„ < 200 keV [FT]. Rele- 
vant calculations for the spectrum of the medium energy 
(£ u ) ~ IMeV neutrinos emitted from thermal processes 
has become available recently 52 , 53 , 54 . Purpose of this 
article is to develop accurate methods and discuss various 
theoretical and practical (important for detection) aspects 
of the neutrino spectra from astrophysical plasma process. 
This could help experimental physicists to discuss possi- 
ble realistic approach to detect astrophysical sources of 
the neutrinos in the future. 

2 Plasmaneutrino spectrum 
2.1 Properties of plasmons 

Emissivity and the spectrum shape from the plasmon de- 
cay is strongly affected by the dispersion relation for trans- 
verse plasmons (massive in-medium photons) and longi- 
tudal plasmons. In contrast to transverse plasmons, with 
vacuum dispersion relation uj(k) = k, longitudal plasmons 
exist only in the plasma. Dispersion relation, by the def- 
inition is a function uo(k) where hu> is the energy of the 
(quasi)particle and hk is the momentum. Issues related to 
particular handling of these functions are discussed clearly 
in the article of Braaten and Segel [TS] . We will repeat here 
the most important features of the plasmons. 

For both types, plasmon energy for momentum k = 
is equal to loq. Value loq = u>(0) is refereed to as plasma 
frequency and can be computed from: 

oo 


where v = p/E, E = \Jp 2 + m\ (H = c = 1 units are 
used), m e ~ 0.511 MeV and fine structure constant is a = 
1/137.036 [55] . Functions /i and j % are the Fermi-Dirac 
distributions for electrons and positrons, respectively: 

= e (E-n)/kT + 1 > h = e {E+n)/kT + l ■ ( 3 ) 

Quantity [i is the electron chemical potential (including 
the rest mass). Other important parameters include first 
relativistic correction iO\ : 

oo 

^ = £/£(§^ + (4) 



maximum longitudal plasmon momentum (energy) k max : 

°° 2 

kLax = u 2 max = — / — ( - ln — — - 1 ) (fi + / 2 ) dp 
ir J E \v 1 + v ) 
o 

(5) 



Table 1. Plasma properties for typical massive star during Si 
burning. All values in MeV. 
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LOA 


0.32 


1.33 


0.074 


0.070 


0.086 0.133 


0.002 



and asymptotic transverse plasmon mass m^: 

oo 

m t=^J P ^{h + h)dp. (6) 
o 

Value m t is often referred to as thermal photon mass. We 
also define parameter u»: 



interpreted as typical velocity of the electrons in the plasma 
|15j . Axial polarization coefficient is: 

oo 

^ = v/é( 1 -i wa ) (/i_/a) *- (8) 



Value of the lua is a measure of the difference between 
neutrino and antincutrino spectra. Set of numerical values 
used to display sample result is presented in Table [JJ 

Values ojo,LJ max ,mt define sub-area of the Lo-k plane 
where dispersion relations for photons LOt(k) and longitu- 
dal plasmons u>i(k) are found: 

max (fc, lo ) < LUi(k) < co max , < k < k max (9a) 

\Jk 2 +loI< Lj t (k) <\Jk 2 +m%, < k < oo (9b) 
Dispersion relations are solution to the equations [15J: 

k 2 =n l (tj l {k),k) (10a) 

k 2 =LU t (k) 2 -n t {LO t (k),k) (10b) 

where longitudal and transverse polarization functions are 
given as an integrals: 

oo 

Aa f p 2 f ui LJi+vk mf -k 2 \ 

o 

(lia) 

OO 

4a f p 2 f lo 2 lo 2 ~k 2 u t . uJt + vk\ 



. ( llb ) 

Typical example of the exact plasmon dispersion re- 
lations (dash-dotted) is presented in Fig. [TJ As solving 
eqns. (jlOai Il0b[) with (fTTj) is computationally intensive, 
three levels of approximation for dispersion relations are 
widely used: 

1. zero-order analytical approximations 

2. first order relativistic corrections 

3. Braaten&Segel approximation 
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Fig. 1. Longitudal and transverse plasmon dispersion relation uji t t(k) for plasma parameters from TableQ] Exact result (dot- 
sahed) is very close to the Braaten & Segel approximation (solid) . Zero-order (dotted) and first order (dashed) approximations 
are very poor, especially for londitudal mode (left). 

Longitudal plasmon In-medium photon 




k [MeV] k [MeV] 

Fig. 2. Longitudal and transverse plasmon mass. Dotted lines on the right panel show asymptotic transverse mass. Line 
dashing the same as in Fig. [T] 



2.1.1 Approximations for longitudal plasmons 



in an implicit form: 



For longitudal plasmons, the simplest zero-order approach 
used in early calculations of Adams et al. [12] and more 
recently in [53j for photoneutrino process is to put simply: 



(12) 



where ujq is the plasma frequency @ . Maximum plasmon 
energy uj max = luq in this approximation. Zero-order ap- 
proximation is valid only for non-relativistic regime, and 
leads to large errors of the total emissivity [12] . 

First relativistic correction to (|12[) has been introduced 
by Beaudet et al. [12]. Dispersion relation u>i(k) is given 



w, 



3 9 ^ 



with maximum plasmon energy equal to: 



3 

—u> 
5 



(13) 



(14) 



This approximation, however, do not introduce really se- 
rious improvement (Figs. [2 [2] (left) Breaking point 

was publication of the Braaten&Segel approximation |15j . 
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Using simple analytical equation: 



3 4 

vi 



U t 



ln 



2v*k ui ~ v *k 



- 1 



(15) 



where v* is defined in ([7]) one is able to get almost exact 
dispersion relation, cf. Figs.Q]&El panels. Solution to 
the eq. (p~5|) exist in the range 1 < k < k^ x , where, in this 
approximation, maximum longitudal plasmon momcntum 
is: 

what gives value slightly different than exact value (Fig.Ø 
left), but required for consistency of the approximation. 



2.1.2 Approximations for transverse plasmons 

For photons in vacuum dispersion relation is ui t = k. Zero 
order approximation for in-medium photons is: 




valid for small k and: 
u 2 = 



k <C ui 



k > ui 



(17a) 



(17b) 



valid for very large k. Formulae (|17a|) and (|17b[) provide 
lower and upper limit for realistic 0J t (k), respectively (cf. Fig 
right panel, dotted). First order relativistic corrections 
lead to the formula: 



Fig. 3. Fenmann diagrams for plasmon decay. 

bisection method. Zero and first-order approximations (|12l 
I17a[ I17b[) with limiting values J9]) provide starting points 
and ranges. Approximation has been tested by [56] and 
is considered as the best available [2D]. Errors for part 
of the kT-fi plane where plasmaneutrino process is not 
dominant may be as large as 5% [56] • At present, these 
inaccuracies are irrelevant for any practical application, 
and Braaten&Segel approximation is recommended for all 
purposes. 



2.2 Plasmon decay rate 

In the Standard Model of electroweak interactions, mas- 
sive in-medium photons and longitudal plasmons may de- 
cay into neutrino-antineutrino pairs: 



I 



7 



(22) 



1 nfc 
5^" 



Ut 



with asymptotic photon mass: 



Finally, Braaten&Segel approximation leads to: 



u. 



hr+wt. 



3u 2 



1 



v 2 k 2 , u t + v*k 
ln 



°2v 2 k 2 \ 2cutv*k Ut — v*k 
Asymptotic photon mass mf s derived from (I20p is 



(-f S ) = 



3l4 
2v 2 



1 - 



1 



2u* 



vi 1 
■ ln ■ 



1 — v t 



(18) 



(19) 



(20) 



(21) 



In the first-order calculations two Feynmann diagrams 
(Fig. [3]) contribute to decay rate [T5"ll5"2"] . 

For the decay of the longitudal plasmon squared ma- 
trix element is: 



Mf 



<~<2 fi2 

^f^v t. ,2 
na v 



2KQ 1 K-Q 2 



K 2 

2 k ■ qi k ■ q 2 

k 2 



Qi-Qa 



(23a) 



where K = (u, k) is four momentum of the plasmon. 
Qi = (£i,qi) and Q2 — (£2,^2) is four-momentum of 
the neutrino and antineutrino, respectively. 

Squared matrix element for decay of the massive pho- 
ton is: 



This is slightly smaller (left panel of Fig. [5] dashed) than 
exact value (solid line). 

All four relations are presented in Fig. [T] Differences 
are clearly visible, but they are much less pronounced 
for transverse than for longitudal plasmons. Inspection of 
Fig.[2]reveals however, that in the large momentum regime 
asymptotic behavior is correct only for exact integral re- 
lations IjlObp and may be easily reproduced using (|17b[) 
with uit from ([6]). 

Let us recapitulate main conclusions. Braaten&Segel 
approximation provide reasonable approximation, as non- 
linear equations (fT5|) and ([20]) are easily solved using e.g. 



Mt = 



Gi 



ira 



(c 2 v n 2 + c 2 A n 2 A ) (e 1 £ 2 - 



k • qi k • q 2 



+ 2C v c A n t n A 



£1 k ■ q 2 - £2 k ■ qi 

k 



(23b) 



where Ilt is defined in (|llb[) and axial polarization func- 
tion II A reads: 



n A = 



2a uji-k 2 



p 2 ( ut_ , Ut + vk 



E 2 \2vk uit — vk u\ — v 2 k 2 



(h- f 2) dp 



(24) 
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Table 2. Relative weight of the M t 2 
neutrinos. 



23b[l terms for e and ji, r 



Flavor 



electron 
mu/tau 



Vector 

n 2 A 



Axial 



Mixed 



For massive photons gt = 2 and for longitudal plasmon 
ffl = l- 

Differential rateio has been derived for the first time 
in [55]. Here, we present result in the form valid for both 



0.74 
0.07 



0.02 
0.39 



2Cv Ca^q^a 

(c v ^+Ca^a ) 2 types of plasmons, ready for calculations using any avail- 
0.24 able form of dispersion relation: 



0.54 



Fermi constant is G F /{hc) 3 = 1.16637(1) x 1CT 5 GeV" 
[55 and, in standard model of electroweak interactions, 
vector and axial coupling constants are: 



d 2 R t 
d£i d£ 2 



4 Z i M if- 

7T 



JiS 



(31) 



C v = -+2sm 2 9 w> C%=- 



2 sin 2 9 



w. 



° A ~~2 



(25) 



(26) 



whcrc i = l or i = t. Product S of the unit step functions 
in (|3Tj) restrict result to the kinematically allowed area: 

S = e(4£ 1 £2-m?)G(e 1 +£2-u )e(u max -£ 1 -£ 2 ) (32) 

Four-momenta in the squared matrix element are: 



Cy Ca do not contribute at all and term C\ x 
pressed relative to the term beginning with C v 



for electron and /i, r neutrinos, respectively. The Weinberg 
angle is sin 2 9 W = 0.23122(15) [55]. 

Terms containing Ca (so-called axial contribution) in 
(|23b[) are frequently treated separately [52] or removed at 
all [3]. In calculations concentrated on the total emissiv- 
ity this is justified as anti-symmetric term multiplied by 

. is sup- 
x ... by 

four orders of magnitude [3]. However, if one attempts to 
compute neutrino energy spectrum all three terms should 
be added together, as mixed V-A ,, channel" alone leads 
to negative emission probability for some neutrino energy 
range (Fig. [6]), what is physically unacceptable. These 
terms remains numerically small but only for electron neu- 
trinos. For /i and r neutrino spectra axial part contributes 
at ~ 1% level due to very small value C V ' T = —0.0376 
while still Ca = —0.5. "Mixed" term leads to significant 
diffcrences between z/ MiT and ^ M;T spectra, cf. Fig. [6] Rel- 
ative contributions of the three transverse " channels" for 
electron and /i, r are presented in Table [5] 

In general, all the terms in the squared matrix element 
(|23b|) should be added. We have only two different spectra: 
longitudal and transverse one. 

Particle production rate from plasma in thermal equi- 
librium is: 



K 

m 2 

cos 9 



k' 



(fi, 0,0, ^) 

(£ 2 , £ 2 sin 8, 0, £ 2 cos 9) 

(£1 + £2 , £2 sin 9, 0, £ 1 + £2 cos ( 

K-K = {£ x + £ 2 ) 2 - k' 2 

k' 2 - f ' 2 _ f ' 2 

2£x£ 2 
^uj- t 1 (£ 1 +£ 2 ) 
= £i + £ 2 



whcre uj. 1 denotes function inverse to the dispersion re- 
lation. Jacobian Ji arising from Dirac delta integration in 
J23 is: 



, _ £i£ 2 duii 

Ja — 



dk 



k=k' 



(33) 



Residue factors Zi are given in ([30]) and (129|) . Maximum 
energy to rnax in (|32[) for longitudal plasmons must be in the 
agreement with particular approximation used for L0i(k): 
ui , ([14]) or ([16]) for zero-order ([12]), first-order ([13]) or 
Braaten&Segel (|15[) approximation, respectively. For trans- 
verse plasmons uj max —> 00 and last function in (|32| has 
no effect and may be omitted. 



Ri 



fj, 



(27T)' 



f /i , 9 d J k d d q 



d 3 k <i 3 qi d 3 q2 ^"^ Longitudal neutrino spectrum 



2£ 2 
(27) 

where i = l for longitudal mode and i — t for transverse 
mode. Bose-Einstein distribution for plasmons / 7 * is: 



f- 



1 



7* 



(28) 



and residue factors Z tt i are expressed by polarization func- 
tions n tA (|llb[[TTa]) : 



2.3.1 Analytical approximation 

We begin with general remark on the spectrum. Note, that 
eq. (I3ip is symmetric for longitudal mode under change 
£ 1,2 — * £ 2,1 because ()23a[) is symmetric with respect to 
exchange Qi$ — * Qa,i- Resulting energy spectrum is thus 
identical for neutrinos and antineutrinos. This is not true 
for transverse plasmons with axial contribution included, 
cf. Sect.EH 



zr 1 = 1 



dilt 
duj 2 

w? m 

k 2 LO 2 



(29) 
(30) 



1 Double differential rate d 2 Ri/d£dcos 8 has an identical 
form as (|31|l but now four momenta cannot be given explicitly, 
unless simple analytical approximation for uji(k) is used. Ana- 
lytical approximations for the specrum shape are derived this 
way. 
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Using zero-order dispersion relation for longitudal plas- 
mons ()12p we are able to express spectrum by the elemen- 
tary functions. Longitudal residue factor Z t is now: 



Z? = 1, 



(34) 



and Jacobian J; resulting from the integration of the Dirac 
delta function is: 



1. 



(35) 



Now, differential rate d 2 Rj dSd cos 9 (cf. (|3Tj) and foot- 
note[T]) becomes much more simple and integral over d cos 9 
can be evaluated analytically. Finally, we get the longitu- 
dal spectrum: 



dR 
~d£ 



G F 2 C v 2 u; 7 f(£/u ) 



1260 tt 4 a fi 3 c 9 e^/ kT 
where normalized spectrum is: 
105 



(36) 



+ 3(1 - 2æ) 2 ln(l - 2x) 



16x 6 



2x z + 6x - 3) 



(37) 



Let us note that / is undefined at x — 1/2; use limit 
instead: 

lim f (x) = 105/32. 

x->l/2 

Function f (x) is symmetric with respect to point x — 1/2, 
where / has a maximum value (Fig. ØJ dotted line). 

In this limit, correct for non-relativistic, non-degenerate 
plasma, average neutrino and antineutrino energy is {£) = 
and maximum v energy is ujq. 

Inspection of Fig. 0] reveals little difference between 
analytical result (f3"6")l and result obtained with first-order 
relativistic corrections to the dispersion relation (fl~3]) . 



2.3.2 Numerical results 

Simple formula (|36|l significantly underestimates flux and 
the maximum neutrino energy, equal to io max rather than 
ojq. Therefore we have used Braaten & Segel approxima- 
tion for longitudal plasmon dispersion relation. 

To derive spectrum we will use form of differential rate 
(|3"Tj) provided by [52]. In the Braaten&Segel approxima- 
tion: 



Z? s = 



2K 



J, BS 



k 2 3u 2 - 
k 2 1 



"f + vi 
Pi 



3u> 2 ( 3lui u>i 
2^1" 



£i£ 2 u>i(3i 
v*k 



\2k 3 



,k k 2 (ujf — v 2 k 2 ) 



2v* 
k 2 



Spectrum is computed as an integral of (|3ip over d£ 2 - 
Example result is presented in Fig. 0] Integration of the 
function in Fig. [5] over neutrino energy gives result in wcll 
agreement with both (30) from [15] and (54) from [52] . 
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Fig. 4. Longitudal plasmon approximate analytical (J36J) neu- 
trino spectrum (dotted), with first-order correction used by 
BPS 12, (dashed), and spectrum computed using |15j disper- 
sion relation (solid). Plasma properties according to TableQ] 



2.4 Transverse plasmon decay spectrum 

2.4.1 Analytical approximation 

Derivation of massive in-medium photon decay spectrum 
closely follows previous subsection. Semi-analytical for- 
mula can be derived for dispersion relations (| 17[) . For dis- 
persion relation (|17b[) transverse residue factor Z t is: 



z- = 1, 

polarization function U t is equal to: 



(38) 



(39) 



and Jacobian resulting from integration of the Dirac delta 
function J t is: 

£l+£2 (40) 



J t 



Approximate spectrum, neglecting diffcrences between 
ncutrinos and antincutrinos, is given by the following in- 
tegral: 



64 ir 4 a h 3 c 9 



P(cos9,£/m t ) dcos9 



2£(l-c 



_1 ex P 

where rational function P(ct,x) is: 

1 + 2(ct - l) 2 (2x 2 - l)x 2 



A/kT 



P = 



x(ct - 1) 2 [1 - 2ct{ct - l)x 2 + 2(ct - l) 2 )x 4 



- 1 
(41) 

(42) 
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where k = 2x + (2x) , x = S/mt, a = mt/kT. Interest- 
ingly, spectrum (|43[) is invariant under transformation: 

£"£ = m 2 /4 




w w max 



£,, [MeV] 



Fig. 5. Transverse plasmaneutrino spectrum computed from 
[15] approximation (solid) with upper (|17b[l and lower (|17a|l 
limits for the dispersion relation (dotted). First-order rela- 
tivistic correction leads to the spectrum shown as dashed line. 
Plasma parameters as in Fig. [4] 



Result presented in Fig. [5j show that spectrum (|41[) 
obtained with dispersion relation (|17b[) agree well in both 
low and high neutrino energy part with spectrum obtained 
from Braaten&Segel approximation for dispersion rela- 
tions. Dispersion relation (|17ap produces much larger er- 
ror, and spectrum nowhere agree with correct result. This 
fact is not a big surprise: as was pointed out by Braaten 
[16 dispersion relation is crucial. Therefore, all previous 
results, including seminal BPS work [12], could be eas- 
ily improved just by the trivial replacement u>o — » m t . 
Moreover, closely related photoneutrino process also has 
been computed [TI l l?l [T71fT4] with simplified dispersion re- 
lation (|17a[) with ujq. One exception is work of Esposito 
et. al. [57j . It remains unclear however, which result is 
better, as accurate dispersion relations have never been 
used within photoneutrino process context. For plasma- 
neutrino, Eq. (|17b[) is much better approximation than 
(I17a|) . especially if one put m t from exact formula ([6"|). 
High energy tail of the spectrum also will be exact in this 
case. 

As formula (l4"Tj) agree perfectly with the tail of the 
spectrum, we may use it to derive very useful analyti- 
cal expression. Leaving only leading terms of the rational 
function (Sl) 



P{ct,x) ~ - c±Y 2 

one is able to compute integral (|4"Tj) analytically: 



\{£) 



Gla 



F^V 



6ATT 4 ah 3 c 9 



k ln | f 

a 



an/2 



1 



(43) 



and all results obtained for high energy tail of the spec- 
trum immediately may be transformed for low-energy ap- 
proximation. The asymptotic behavior of (|4"3"| for £ ^> kT 
is of main intcrest: 



\{£ v ) =AkTm 6 t exp 
where for electron neutrinos : 

GpCy 1 



kT 



(44) 



A = 



87r 4 a h A é> 



= 2.115 x 10 JU [MeV^cm^s -1 ] 



and m t , kT are in MeV. For \i, r neutrinos just replace A 
with A(C$ T /Cp) 2 . 

Formula (|44[) gives also quite reasonable estimates of 
the total emissivity Qt and mean neutrino energies {£ v ): 



Qt = A kT 3 m\ 
(£u) = kT 



(45a) 
(45b) 

For a comparison, Braaten & Segel [15] derived exact for- 
mulae in the high temperature limit kT ^> ujq: 



Qf 



BS 



g|qR(3) 

127r 4 a 



kT 3 m\ = 0.8 A kT 3 m\ (46a) 



(S^ s ) = ^-kT = 0.73 kT 



(46b) 



Formulae above agree with ~25% error in the leading co- 
efficients. 



2.4.2 Numerical results 

Calculation of the spectrum in the framework of Braaten&Segel 
approximation requires residue factor, polarization func- 
tion [K] (transverse&axial) and Jacobian [52] : 



Z. 



BS 



n 



BS 



3 W g W ? + ( W ? + fe2)( W ? - vi fc2) _ 2uH^ ~ fe 2 ) ' 

(47) 
(48) 



3wg 



vik 2 



u t . u t + v*k 
ln 



n 



BS 



w a k 



k 2 2w*fe uit — v*k 
k 2 3LU 2 -2{uj 2 ~k 2 



vik 2 



T BS _ 



oBS _ 9a; 

Pt Avlk 2 



1 



£i£ 2 

1 (v*k 
6 l w ( 



BS 



1 



W Pt 



iss 



3w* \ ut 

v*k I ut — v*k 



(49) 
(50) 
(51) 



Example spectrum, computed as an integral of (|3ip 
over d£ 2 is shown in Fig. [5j 
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t; 4 



y 3 



co 

^2 
CC 



w w max 



0.5 



£,, [MeV] 



Fig. 6. Spectrum of the muon neutrinos (dotted) and antineu- 
trinos (dashed) from transverse plasmon decay. Contributions 
to the spectra from so-called mixed „vector-axial channel" pro- 
duces significant differences. For electron flavor, contribution 
from "mixed channel" lead to unimportant differences. For 
both flavors contribution from "axial channel" remains rela- 
tively small: 10 -4 for u e and 10~ 2 for v^. Overall contribution 
to the total emissivity from /i, r flavors is suppressed relatively 
to electron flavor by a factor (C^ T /Cyf ~ 3.3 x 10" 3 . 



3 Summary 

Main new results presented in the article are analytical 
formulae for ncutrino spectra (|361 141[) and exact analyti- 
cal formula (|44|) for the high energy tail of the transverse 
spectrum. The latter is of main interest from the detection 
of astrophysical sources point of view: recently available 
detection techniques are unable to detect keV plasmaneu- 
trinos emitted with typical energies (£„) ~ w /2 (Fig. [U 
[5]), where ojq is the plasma frequency ([2]). Tail behavior 
of the transverse spectrum quickly "decouple" from loq 
dominated maximum area, and becomes dominated by 
temperature-dependent term exp (— £ v /kT). Calculation 
of the events in the detector is then straightforward, as de- 
tector threshold in the realistic experiment will be above 
maximum area. This approach is much more reliable com- 
pared to the typical practice, where an aver age neutrino 
energy is used as a parameter in an arbitrary analytical 
formula. 

Analytical formulae for the spectrum are shown to be 
a poor approximation of the realistic situation, especially 
for longitudal plasmons (Fig.ØJ). This is in the agreement 
with general remarks on the dispersion relations presented 
by Braaten [16] . On the contrary, Braaten & Segel [15] ap- 
proximation is shown to be a very good approach not only 
for the total emissivities, but also for the spectrum. Excep- 
tion is the tail of the massive photon decay neutrino spec- 



10 2 



1(r 

J 10 2 
> 

CD 



<o 

I 10 18 



10 1 



10" 



10" 



1 <>><S>>max 

£„ [MeV] 




Fig. 7. Typical spectra from the plasma process. Dotted 
line is a longitudal and dashed transverse spectrum. Only 
~ exp(— £„/A;T) tail of the transverse spectrum (solid line) 
contributes to (possibly) detectable signal. Plasma properties 
according to TableQ] 



trum: Braaten & Segel [15] formulae lead to underestimate 
of the thermal photon mass while the formula (|44[) gives 
exact result. Numerical difference between mt from ([6]) 
and pT]) is however small [TS] . Calculating of the emissiv- 
ities by the spectrum integration seems much longer route 
compared to typical methods, but we are given much more 
insight into process details. For example, we obtain exact 
formula for the tail for free this way. Interesting surprise 
revealed in the course of our calculations is importance 
of the high-momentum behavior of the massive photon. 
While mathematically identical to simplest approach used 
in the early calculations, formula (|17bj) gives much better 
approximation for the total emissivity than ()17a|) . 
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